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ABSTRACT 


The analysis of sound propagating by multiple paths in 
an ocean at short ranges has been conducted using a Modified 
Time Delay Spectrometry (TDS) technique. In this version of 
TDS, a source driven by a linear FM slide and an HP3561A 
Dynamic Signal Analyzer are used to measure the amplitude as 
a function of frequency of signals traveling by different 
paths and having different arrival times. Two sets of data 
from the acoustic test ranges at the Naval Undersea Weapons 
Engineering Station were analyzed for different environmen- 
tal conditions to determine the relative amplitudes of the 
directly propagating and surface reflected signals. Compar- 
isons with simple rough surface scattering theory showed 
reasonable agreement. Results and control software are 
presented and discussed. Recommendations for future 


applications are made. 
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The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases 
of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification is at the eilsk or seneeuscer- 


ae. 


i i. 


ey. 


INTRODUCTION ---------------- - - - - - - - - - - - - - - - - 
A. BACKGROUND ---------------------------------- 
B. OBJECTIVES ----------------------- -- -- - ------ 
Peeper ockherT FON AND CONS IDERATIONS ---------- 
A. THE MULTIPATH PROBLEM ------~----------------- 
1. Assumptions ----------------------------- 
2. Multipath Geometry ---------------------- 
3. Surface Reflection ---------------------- 
4. Environmental Considerations ------------ 
B. TDS CONSIDERATIONS -------------------------- 
ine The LFM Pulse --------------------------- 

2. Response of the HP3561A Dynamic 
Signal Analyzer ------------------------- 
TDS MULTIPATH MEASUREMENT SYSTEM ---------------- 
A. RANGE AND TEST CONFIGURATION ---------------- 
B. THE DATA PROCESSING SYSTEM ------------------ 
1. Hardware -------------------------------- 
2. Software -------------------------------- 
a. Quick Analysis ---------------------- 
b. Normal Analysis ----------<-<---------- 
fae eee 
ad. Analysis Products ------------------- 
DATA ANALYSIS ------------ 9-2-2 
A. ANALYSIS PARAMETERS ------------------------- 


TABLE OF CONTENTS 


10 


10 


iE 


LZ, 


iz 


lea 


ie 


18 


Zu 


Z2 


Ze, 


24 


Ze 


Ze) 


28 


28 


20 


a2 


5 


36 


26 


ey 


oy 


1. Time Delays 


2. The Rayleigh Reflection Coefficient ----- 


3. Selecting Analyzer Settings ------------- 


B. MEASUREMENTS 


C. VARIABILITY OF RESULTS -------92-92----------- 


Ve DISCUSSION AND CONCLUSIONS ---------------------- 


A. DISCUSSION 


ile Procedure 


ee Software 


B. CONCLUS#eNs 


APPENDIX A: RANGE AND TEST DATA ---------------------- 
APPENDIX B: TABULATED MEASUREMENT RESULTS ------------ 
LIST OF REFERENCES -------2 3-92-2222 


INITIAL DISTRIBUTION 


Biot 


ay 


oe 


40 


42 


44 


47 


47 


47] 


49 


50 


BZ 


a3 


a9 


60 


PESO r TABLES 


DISTANCES AND TIME DELAYS ----------3---- rrr err 


Peon OoenONtwotNol lt vVittihs AND ENSERTION LOSSES <--=<=- 


Peover tie eCOn ey lerr hl At Sse Leer PREOUENCIES ----- 


oie rottCAab COMPARTESON OF SURFACE REFLECTED 


LEVELS 


39 


43 


44 


46 


Pol OF TReURES 


Simple Multipath Geometry ------------------------ 
The LFM Pulse ------------------------------------ 
TDS Hardware ------------------ - - - - 
Structure of the TDS 1 Program -<-<<<<-<<<<<<<=———-— 


Storage of an LFM Pulse in the HP3561A Dynamic 
Signal Analyzer Time Buffer for a Typical Case --- 


Sound Velocity Profile (Nanoose-February) -------- 
Quick Analysis 3-D Spectral Map ------------------ 


Quick Analysis and Normal Analysis Products -----= 


14 


Zo 


Ze 


mk 


a0 


ac 


41 


45 


Pa ACKNOWLEDGEMENT 


Partial support and the data necessary for the 
successful completion of this thesis from the Acoustics 
Division of the Research and Engineering Department at 
NUWES, Keyport, is gratefully acknowledged. I want to thank 
Mr. John Burwell for the dedicated support and time he has 
given to this project. I also wish to thank the faculty and 
staff at NPS whose instruction has had a significant effect 
on shaping my understanding of Underwater Acoustics. In 
particular, a special thanks to my thesis advisor, Prof O. 
B. Wilson, for his guidance in this endeavor. 

Finally, to my most staunch supporter whose love 


extended to helping me type my thesis, my wife Lynette. 
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Ps BACKGROUND 

The Naval Undersea Warfare Engineering Station (NUWES) 
in Keyport, Washington, performs measurements on torpedo 
radiated noise as part of the test and proofing process of 
torpedoes. From these measurements torpedo radiated noise 
baselines are established, torpedoes are accepted or 
rejected on the basis of meeting established radiated noise 
criteria, and design improvements in noise reduction can be 
evaluated. For these reasons, it is important that sound 
transmission loss between the source, the torpedo, and the 
receiver, a system called the Noise Recording System (NRS), 
be accurately modeled. The current model, based on 
empirical data, appears to be satisfactory under normal 
conditions. However, under unusually rough surface 
conditions scattering effects may disturb the measurement. 

A limited number of experiments and theoretical studies 
have been conducted by Naval Postgraduate School (NPS) 
Students and Faculty to determine the significance of rough 
surface reflection in the multipath acoustic transmission 
loss problem. Brekke [{Ref. 1] developed a computer- 
controlled FFT-based dynamic signal analyzer variation of a 
technique called Time Delay Spectrometry (TDS) to measure 


the separate contributions of multipath propagation. 
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TDS, originally developed by Richard C. Heyser in 1967 
(Ref. 2], utilizes a transmitted Linear Frequency Modulated 
(LFM) pulse or "swept tone" with constant amplitude. A 
frequency tracking spectrum analyzer uses the difference in 
arrival times to discriminate between the acoustic signal 
traveling by a direct path and those signals traveling by 
reflected paths. 

In a thesis at NPS, LT Ward {Ref. 3] outlined the multi- 
path problem including environmental considerations for the 
NUWES test ranges at Dabob Bay and Nanoose. Further, he 
generated a computer model for predicting transmission loss 
under varying sea states. This model has not yet been veri- 
fied experimentally. However, TDS may now make this 


possible. 


ba OBJECTIVES 

The objective of this thesis is to apply the modified 
TDS multipath measurement system to measure the effect of 
surface roughness on propagation loss as a function of 
frequency and compare the results with theory. The problen, 
theory, and special considerations will be presented first, 
followed by a description of the TDS systems. Analysis of 
data taken at the test range at Nanoose will then be 
presented including a brief discussion of the variability of 
the results. Results, conclusions and suggestions for 


future research will also be discussed. 
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Il. PROBLEM DESCRIEPTIONSANDSCCNS TDERAEeiiS 


A. THE MULTIPATH PROBLEM 

For the measurement of radiated noise, during part of 
the test run, the torpedo passes by a vertical array of 
three calibrated omnidirectional hydrophones. Maximum 
horizontal range between the torpedo and array is less than 
1000 yards. Torpedo and hydrophone arrays are at 
approximately mid water depths. Water depth is typically 
about 600 feet in Dabob Bay and about 1200 feet at the 
Nanoose range. Transmission loss due to absorption can be 
neglected since geometrical spreading losses are between 40 
and 60 dB greater than absorption at these ranges at sound 
frequencies of interest between 50 Hz and 30 kHz. 

Since the tests are conducted in relatively shallow 
water, acoustic reflection from the boundaries can make 
Significant contributions to the total acoustic intensity at 
the receiver. Therefore, in order to model the multipath 
problem adequately, geometry, boundary conditions, and 
environmental effects must be considered. 

1. Assumptions 

It is assumed that the acoustic signal is of small 
amplitude and that sound propagation can be described by the 
linear inhomogeneous wave equation. Thus the water channel 


acts as a linear filter. It is also assumed that the water 


ae 


column containing the source and the receiver is 
homogeneous. The bottom is considered to be flat and 
iemuzontal, 

2. Multipath Geometry 

Figure 2.1 depicts simple multipath geometry. 
According to the approach developed in Albers [Ref. 4:pp. 
49-51] and presented by Ward [Ref. 3:pp. 87-88] and Brekke 
(Ref. l:pp. 25-27], the acoustic signal at the receiver is 
the sum of source signals traveling the direct path and the 
contributions of "image" sources traveling various reflected 
paths. 

Following Brekke's notation [Ref. l:p. 27], the 
source depth is ZS and the receiver depth is ZR. The 
horizontal separation between source and receiver is R and 
the channel depth is H. In terms of these quantities, the 
direct path distance (XR), the surface reflected path 


distance (XSR) and bottom reflected path distance (XBR) are 


given by: 
XR = [R2 + (ZR - 2S)2]1/2 (2-1) 
XSR =[(R2 + (ZR + ZS)2]1/2 , (2-2) 
and 
XBR = [R2 + (2H - ZR - 2S)2]}1/2 . (2-3) 
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Figure 2.1 Simple Multipath Geometry 
(Extracted from Ref. 1] 
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For periodic signals, the contributions of the image 
sources must be considered with respect to both relative 
phase and amplitude. The phases of the reflected signals 
are shifted with respect to the directly propagating signal 
by an amount proportional to the difference in propagation 
path. Additionally the surface reflected signal undergoes a 
complete phase reversal at reflection. The delayed arrival 


times, for a given sound velocity c, are computed as 


follows: 
Surface time delay is tac= ee (2-4) 
Bottom time delay is th = aR (22s) 


and the time delay between arrivals of surface and bottom 


Pome eccdmcounds, OF Intcer—-reflected time delay t,, is 


~ XSR_- XBR (2-6) 
r C 


ie 


Assuming an omnidirectional source radiating simple 
harmonic waves, the acoustic pressure wave function at a 
range r is given by: 


P = (A/r) exp j(wt-kr) (27 ) 


where A is the pressure amplitude at unit distance, k is the 


wave number W/c), and is the angular frequency [Ref. 5: 


LS 


p. 112]. The acoustic pressure at the receiver due to the 


Signal arriving via the direct path is given by: 


Pq = (A/XR) exp j(wt-(K{XR}) ] (2-8) 


The contribution of the surface "image" is given by: 


P, = -R, (A/XSR) exp j[uwt-(k{XSR}) ] (2-9) 


where Rs is the coefficient of reflection for the surrcee 
(Ref. 3:p. 90]. Similarly, the contribution due to bottom 


reflection 1S given by: 


Ph = Rp (A/XSR) exp) [wt-(k{XBR}) ] (2-10) 


where R, is the coefficient of reflection for the bottom: 

It is possible that Ry, may have an additional phase term 
determined by the nature of the water-bottom composition and 
the angle of incidence (Ref. 3:p. 90]. Assuming that the 
direct and reflected waves are coherent one finds the total 
acoustic pressure at the receiver is the sum of several 


CONCEIDULIONS: 


Ptot = Pq +t Pg + Ph - (2-105) 


Equation (2-11) can be rewritten as: 


Prot = M Pq (2-12) 
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where 


ies exp (jwt.) +yes R, exp(jut,)] ‘ (2-13) 


Maisie XSR 6S 


Merdmbret.e3:p. 948) defines this factor as the "Multipath 
effect." 

If the phases of the signals are incoherent, the 
intensity of the combined signals is the sum of the 


intensities of the separate signals, and is given by: 
I+ot = Ig t+ Ig t1Ip : e214 ) 
Equation (2-14) can then be rewritten as: 
ITtot = Mr Iq _ , (2=15) 
where 
My = (11/7 tae 1,/ 14a) : (2-16) 
Equation (2-16) is defined as the Intensity Multipath 
coefficient. The surface reflection intensity, I,, is 
determined by the condition of the surface. 
The reflection coefficients, R, and Rp, are 
determined by the physical properties of the boundaries. Rp, 


is small at both ranges [{Refs. 1, 3] and should be 


IL 7 


considered significant only when either source or receiver 
is located near the bottom. For purposes of this research 
Rp is assumed to be negligible. The surface is assumed to 
provide no mechanism for absorbing acoustic energy and, 
since the specific acoustic impedance of air is small 
compared to that of seawater [Ref 5:p. 127], the energy in 
an acoustic signal incident at the surface is assumed to be 
perfectly reflected or scattered back into the water. 
3. Surface Reflection 

Urick [Ref. 6:p. 128] describes the sea surface as 
"both a reflector and a scatterer of sound and has a 
profound effect on propagation . . . where source and 
receiver lie at shallow depth." Surface "roughness," 
defined in terms of signal wavelength, determines the nature 
of the acoustic signal after interaction with the surface. 
Urick defines the Rayleigh parameter, KH sin 6, as the 
criterion for acoustic "roughness," where k is the wave 
number (w/c), H is the rms wave height (peak to trough) and 
is the grazing angle [Ref. 6:p. 129}. When kH Sin 6 >> 1, 
the reflection process is called scattering and the sea 
surface randomly scatters all the energy from an incident 
acoustic signal. For kH sin 96 << 1, the surface behaves as 
an acoustic mirror and the reflected wave is completely 
coherent with the incident wave [Refs. 3, 5, 6]. 

The vertical elevation of the sea surface is a time- 


variant random process and is usually modeled statistically. 
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A relatively good first approximation of the probability 
distribution of the surface wave height [Ref. 3:p. 66], 
particularly when kH sin 6 < 1 [Ref. 7:p. 344], is the 
stationary Gaussian process. According to developments of 
Beckmann and Spizzichino [Ref. 8:pp. 73-74], the mean time 
average scattering coefficient <q> for a normally 


distributed surface 1s given by: 


<q> = X(Kz) do (2-17) 


where q, is unity (smooth surface), X(K,) is the spatial 
Fourier transform of the surface probability density 
Function or "characteristic function," and kz is the wave 
propagation vector for the acoustic signal. The 


characteristic function for a Gaussian distribution is given 


Dy: 

X(kz) = exp (- of kz,%/2 ) (2-18) 
where o 1s the mean-to-peak rms wave height. The 
propagation vector, k,, for the case where the reflected 


angle ¢%, is equal to the incident angle 4%;, is given by: 


kz, = 2k cos 64 3 (2-49) 


io 


The term kz, 1s equivalent to the Rayleigh parameter. 
Beckmann and Spizzichino define it as g!/2 and use it as a 
measure of surface roughness [Ref. 8:p. 82]. Using the 
definition of g and Equation (2-15), the mean time average 


scattering coefficient of Equation (2-13) becomes: 


<q> =sexp(-G72) ; (2 
The power reflection coefficient, 15 is the ratio 


of acoustic energy scattered by a rough surface to the 
energy scattered by a smooth surface (unity) [Ref. 3:p. 68] 


and is defined by the expression: 


R = <qq*> = exp (-4k2 o* cos? O;) . (2-2 15) 
T 
The surface amplitude reflection coefficient, R,, is then 


Givensby [Ref sesspamees |]: 
ee Bee = exp(-2k* g@cos* 4) . (2-22) 


Based on experimental determinations, Clay and Medwin [Ref. 
7:p. 344], show that for values of gi/2/2 greater than 
unity, Equation (2-19) tends to under-estimate the true 
value of Rg. They cite shadowing effects and variation in 


actual wave shapes as reasons for these differences. 
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If it assumed that there is negligible bottom 


reflection, the multipath term of equation (2-12) becomes: 
a 28h mes = 
Ms = [1 exp ( jut.) ] : (2-23) 


The sound pressure level at the receiver due to sound 


travelling by the direct path is given by: 


SPLp = SPLy - 10 logzo |Me!? , (2-24) 


where 


SPLy = 20 logyo |Ptot! - (2-25) 
Equation (2-22) is the total sound pressure level at the 
receiver. 
4. Environmental Considerations 

Brekke [Ref. 1] and Ward [Ref. 3] include extensive 
environmental data for both NUWES test ranges. The key 
environmental variable affecting multipath propagation is 
extreme variation in the sound velocity profile (SVP). 
There is the possibility of strong negative gradients 
occurring above 100 ft in depth. The resulting refraction 
may reduce delay time between the signal traveling the 


direct path and reflected paths so that they may not be 


21 


discernible [{Ref. l:pp. 27-37] or may affect the actual 


angle of incidence at the surface. 


B. TDS CONSIDERATIONS 

An analytical discussion of the modified TDS technique 
1s presented by Brekke [Ref. l:pp. 16-24]. The points of 
the verification to be considered here are the linear sweep 
or LFM pulse and the response of the HP3561A Dynamic Signal 
Analyzer. 

1. The LFM Pulse 

Figure 2.2 represents the LFM pulse in the time 

domain. The amplitude of the signal is constant and for 
convenience is assigned a value of unity. The frequency 
range F is the difference between the start frequency fj and 
the stop frequency f5. The carrier frequency f. 1s defined 
as the average of f5 and ff}. The sweep has a period of T 
seconds before it repeats. In terms of these parameters, 
the instantaneous voltage of the LFM pulse can be determined 


as a function of time in the expression [Ref. l:p. 17]: 


X(t) = cos({(n(F/T)t? - TFt + 2 TEQt] . (2-26) 


The phase of the Signal is: 


o(t) = «a(F/T)t? + 2) Sp ieee (2-27) 
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Figure 2.2: The LFM Pulse 
CEseractea from kerf... 1} 


The instantaneous frequency f; is obtained by taking the 
derivative with respect to time of Equation (2-27) and 
feraing by 27 : 


fj = (F/T)t + fy . (2-28) 


The sweep rate (S) is defined as the derivative with respect 


to time of Equation (2-28): 
S= F/T. (2-29) 


In the context of multipath propagation, the 


instantaneous frequency of the signal propagating by the 
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direct path differs in frequency from the reflected signals 
by an amount equal to the sweep rate, S, times the delayed 
arrival time of the reflected signal given in Equations (2- 
4) angimg2—oon 

2. Response of the HP3561A Dynamic Signal Analyzer 

In the current experimental configuration, the LFM 

pulse is applied to a broad band-width projector. The 
acoustic wave propagates by one or more paths through the 
water to a hydrophone. The electrical signal from the 
hydrophone is recorded in analog form on magnetic tape. The 
recorded signal is subsequently analyzed using the HP3561 
Dynamic Signal Analyzer. The meaSurements are based on a 
1024 point FFT algorithm. The displayed spectrum contains 
400 bins. Frequency resolution of the analyzer is 
determined by the selected frequency span F of the analyzer 


in terms of the bin separation b: 

b = F/400 Hz. (2-30) 
For a frequency span of 5 kHz, the corresponding frequency 
resolution is 12.5 Hz. The time record length, the 
observation time t,p. during which the analyzer receives 


data from the signal, is given by: 


ieee Io “eee. (2-31) 
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This is the minimum amount of time required by the analyzer 
tO obtain sufficient samples to carry out the FFT algorithm. 
The energy contained in each of the 400 bins is calculated 
and made available for display by the analyzer during the 
observation time. In most operating modes, the HP3561A 
displays only the most recent time record data. In Time 
Capture mode, however, the analyzer acquires and stores data 
contained in 40 contiguous time records. Time Capture will 
be discussed in more detail in the next chapter. 

The consequence of these relations is that the delayed 
arrival times given by Equations (2-4) and (2-5) must be 
greater than tope in order to resolve the directly 
propagating and reflected signals [Ref. l:p. 48]. For delay 
arrival times less than tops, the difference between the 
instantaneous frequencies of the direct and reflected 
Signals is less than the resolution of the analyzer and 
thus, the energy of both direct and reflected signals are 
contained in the same bins. 

The possibility of energy leaking into an adjacent 
bin and creating additional measurement errors is reduced by 
using a "window function" feature on the HP3561A. The 
result of using a window function is a minimum bandwidth 


separation, f band, given by [Ref. l:p. 48]: 


f£ band = xr ,; (2232) 


Zo 


where C, a constant determined by the window function 
selected, is expressed in terms of the frequency span. For 
a Hanning window, C is 0.00375. Thus, if the selected 
frequency is 5 kHz, the resulting £ band is 18°75 Hz > =ae, 
in order to resolve the direct and reflected signals, it is 
necessary to insure that the delayed arrival times are 
sufficient to maintain the minimum bandwidth separation. 
This is accomplished by setting a lower bound on the sweep 


rate, S,, given by (Reb. 1:p. 92om. 


om > £ Daneyer delay =, (2-335) 


where t delay is given by Equations (2-4) through (2-6). 


For surface reflection, t_delay = tg. 
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tee 2b Mun! ) PATH OMEASUREMENT SYSTEM 


The TDS Multipath Measurement System 1S comprised of two 

major subsystems [Ref. 1l:p. 39]: 

(1) Range and test configuration 

(2) Data processing system. 
The range and test configuration described here is specific 
to the data used in this research. The test format was 
developed by John G. Burwell of the Acoustics Division, 
NUWES, Keyport, and has been implemented by Keyport 
personnel in previous TDS research. The data processing 
system presented by Brekke [Ref. l:pp. 42-61] has been 
modified. A description of the hardware configuration is 
provided here. The software used is a Slightly modified 
version of the TDS software documented by Brekke [Ref. 1]. 
The key features of the TDS programs will be discussed here 


Bmcescotktware modifications will be documented. 


A. RANGE AND TEST CONFIGURATION 

The data used in this study were obtained at the NUWES 
range at Nanoose, Canada, in February 1986. The acoustic 
projector for the test was the Sonar Acoustic Test System 
(SATS). The receiver was the Noise Recording System (NRS) 
with three calibrated hydrophones at depths of 330 ft., 390 


ft., and 430 ft. Horizontal separation between SATS and NRS 


Ze 


was approximately 300 yds. Wave heights were estimated at 1 
to 2 ft peak-to-peak at the time of the measurements. 

The test signal was a linear sweep from 50 Hz to 30 kHz. 
Three tests were conducted with sweep periods of 1 second, 2 
seconds and 4 seconds. Thirty sweeps were transmitted for 
each test. 

The sound signals arriving at each of the NRS 
hydrophones were recorded on separate channels of a magnetic 


tape recorder. 


B. THE DATA PROCESS@ENG SYS@EM 
1. Hardware 

The equipment configuration used is shown in Figure 
3.la. It is essentially the same as that used by Brekke 
(Ref. l:p. 42-43], with an HP thinkjet printer replacing the 
HP2671G printer. For the routine tests, the graphics and 
printing requirements could be satisfied using the HP 
thinkjet as shown in Figure 3.1b. 

The measurements of the signal are carried out by 
the HP3561A Dynamic Signal Analyzer on the recorded signal 
under computer control. Previously, triggering the Dynamic 
Signal Analyzer was done using a 1 pulse-per-second signal 
derived from the output of the time code generator [Ref. 1}. 
This was not possible, since the period of the LFM pulses 
for two of the tests exceed 1 second. To accommodate the 


longer sweep periods, the trigger signal of an HP3314A 
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Figure 3.1 TDS Hardware 
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Arbitrary Function Generator was set to 4 seconds and 
recorded on the edge track of the data tape. 
2. sottware 

Brekke [Ref. 1] provides a detailed documentation of 
the TDS 1 program. Here, a general discussion of the 
capabilities of the TDS software will be presented with 
emphasis placed on the features used in this research. The 
modifications used in the current version 3C of the TDS 1 
program will also be documented. 

The original TDS program, TDS _ 1, was designed to 
control the HP3561A Dynamic Signal Analyzer in performing 
TDS Multipath Measurements. It 1S written in HP Basic 
program language 2.0 with extension AP2 1. The program is 
structured in a modular design illustrated in Figure 3.2. 
It was designed to analyze a series of LFM pulses sweeping 
from 500 HZ to 20 kHz with a repetition period of 1 second. 
The Set analyzer subprogram in version _ 1 set the HP3561A to 
measure an LFM pulse corresponding to these parameters. In 
order to accommodate the sweep rates and sweep periods for 
the new data described in part A of this chapter, the 
program required some small modifications. The set analyzer 
subprogram which services all the major subprograms was 
altered to allow an operator interactively to select the 
frequency span, start frequency and trigger method of the 


analyzer prior to starting any analysis. 
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The signal is analyzed by one of three methods: 
Quick Analysis, Normal Analysis, and Auto Analysis. 
Referring to Figure 3.2, one sees that these three methods 
form the three major subprograms of TDS 1. These analysis 
methods are discussed in the following paragraphs. 

a. Quick Analysis 

In this method the HP3561A is operated in the 
"Time Capture Mode". The analysis is based on data taken 
from 40 time records stored in the time buffer of the 
analyzer memory. Each time record consists of the data 
processed during one observation time tpp.q-. The total time 
during which the analyzer processes data is 40 x ee 
Figure 3.3 illustrates the storage of the LFM pulse in the 
time buffer. The data are displayed as a 3-dimensional 
magnitude map. Each map consists of 60 traces. The overlap 
(or increment) processing selected determines the time 
increment represented by each subsequent trace. For the 20% 
overlap processing used in the TDS program, each trace 
represents a 0.2 X tope increment of the delay time. The 
analyzer uses 20% new data taken from the subsequent time 
record for each new trace. 

The purpose of the overlap processing is to 
optimize the useful information displayed by the analyzer. 
In selecting an appropriate overlap, the sweep rate, the 
analyzer frequency span, and the frequency span of the LFM 


pulse must be considered. For example, the 4 second LFM 
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pulse sweeps through a 5-kHz frequency span in approximately 


678 ms. The time delay between the direct and reflected 


Signal is about 100 ms. For the 20% overlap processing used 
in the TDS program each trace represents 0.2 x 80 ms = 16 


ms. Thus, the analyzer uses about 50 traces to display 
direct and reflected signals. If, however, 30% overlap 
processing is used, each trace then represents 24 ms. [In 
this case only 35 traces are used to display direct and 
reflected signals, leaving 25 traces displaying no useful 
data. If 10% overlap processing is selected, each trace 
represents 8 ms. The number of traces required to display 
direct and reflected signals is approximately 90 which 
exceeds the 60 traces that can be displayed by the analyzer. 

Quick Analysis provides the operator a quick 
overview of a "Spectral Segment" of a single sweep of the 
Signal to be analyzed. A "spectral segment" is the portion 
of the LFM pulse that can be displayed by the analyzer. If 
the start frequency of the analyzer is 500 Hz and the 
frequency span selected is 5 kHz, then only the "segment" of 
the LFM pulse between these two frequencies can be 
processed. 

From the spectral 3-D map, the arrival times of 
the "Spectral Segment" for both direct and reflected signals 
can be determined by adjusting the start time of the 
analyzer. The time when the reflected signal vanishes from 


the "Spectral Segment" should also be determined. These 
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values are used aS input parameters for Normal Analysis and 
Auto Analysis. After the determination of the arrival times 
and “upper limits of the delayed start time", or time when 
the reflected Signal vanishes, the program displays 
individually all magnitude maps between these time limits 
for operator input. The operator selects peaks for direct 
and reflected signal using the cursor on the Spectrum 
Analyzer. These values are recorded by the program and are 
printed out at completion of the analysis. 

It should be noted here that Quick Analysis 
requires only a single trigger since only one sweep is 
captured. The requirement for a sequence of synchronous 
triggers exists when Normal Analysis or Auto Analysis will 
be performed. 

b. Normal Analysis 

In this method 8 sweeps are RMS averaged to 
reduce effects of signal fluctuation and increase signal to 
noise ratio before making measurements. The measurements 
are made in a manner similar to those in Quick Analysis in 
which the operator must select peak values using the cursor 
on the HP3561A. The program increments the delay step after 
each measurement to measure the frequency bin of the later 
arrival time. These measurements are repeated for every 
delay step in the time interval determined by the 


Quick analysis method. 
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The delay step should increment each measurement 
by a frequency larger than the frequency resolution of the 
analyzer. 

c. Auto analysis 

This method automates the measurements performed 
by Normal Analysis method. Due to a high noise level at the 
lower end of the spectrum, the operator's observation was 
necessary to discriminate between direct and reflected 
Signals. Consequently, the Auto Analysis method was not 
used in this study. 

dad. Analysis Products 

Each of the methods produces a table of voltage 
levels for direct and reflected signals as well as a noise 
level. The program enables the operator to plot the 
magnitude of direct and reflected signals as a function of 
frequency. This requires that all the tabulated frequency 
values be monotonically increasing. This is not always the 
case, since often the reflected path peak is recorded at a 
lower frequency than the previous measurement due to the 
uneven nature of boundary reflection. Using a longer delay 
time will reduce this problem, but with the result that 


fewer samples are recorded. 
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IV. DATA ANALYSIS 


The analysis of a "spectral segment" of LFM pulse with 
the 4 second period is presented here. Procedures used here 
are the same for the analysis of any "Spectral segment" of 
any LFM pulse. Multipath problem considerations are 
discussed with respect to TDS parameters. Measurements made 
using the Normal Analysis method are used to determine the 
effects that surface reflection or scattering have on the 
acoustic field at the receiver. Measurements by the 
Quick Analysis method on several different sweeps are 
compared to the Normal Analysis results to estimate the 


Sweep to sweep variability in results. 


A. ANALYSIS PARAMETERS 
1. Time Delays 

The sound velocity profile given in Figure 4.1 is 
typical for the Nanoose range in February [{Ref. 9]. Because 
the effects of refraction due to this profile are small for 
horizontal ranges of a few hundred yards, the water column 
will be assumed to have a constant sound speed of c = 4850 
ft/sec. The multipath parameters for the problem are given 
in Table 1. 

The minimum time delay given by Table 1 is 82 ms. 
The minimum frequency span setting for the HP3561A is 5 kHz 


corresponding to type = 80ms. 
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Figure 4.1: Sound VETOecrry Profle (Nanoose-February) 
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TABLE 1 


DISTANCES AND TIME DELAYS 


RCVR (depth) Me Cycle) XSR (yds) t. (ms) 
upper (330) 100 yds 233 yds 82 ms 
middle (390) 104 yds 251 yds 91 ms 
lower (430) 210 yds 266 yds 96.4 ms 


Seumce depth: 300 4fte 
Horizontal Range: 100 yds (300 ft) 


c = 4850 ft/sec. 


2. The Rayleigh Reflection Coefficient 


Assuming transmission loss to be due to. only 
~spherical divergence, the sound pressure level at _ the 
receiver due to the direct signal is given by: 


SPLg = SPL{P(1)] - 20 log XR . (4-1) 


The sound pressure level due to a surface reflected signal 


is assumed to be given by: 


SPLp = SPL [P(1)] - 20 log XSR + 10 log (Rg) . (4-2) 


P(1) is the sound pressure at unit distance from the source. 


pubtracting (4-2) from (4-1) yields: 


a 


10 log Rg = SPLp-SPLp + 20 log (XSR/XR) . (4-3) 


If all terms on the right hand side of (4-3) are 


known, then a value for R, can be determined. Using 
equation (2-21), one can solve for the mean-to-peak wave 
height. For six sweeps analyzed using the Quick analysis 


method, these calculations gave an average value for o of 
0.54 ft at frequencies less than 1200 Hz. The standard 
deviation at these frequencies was on the order of 0.05 ft. 
At frequencies above 1200 Hz, the calculated values for O 
decreased. 

These results are consistent with data given by Clay 
and Medwin (Ref. 7:p. 344]. Since the wave heights were 
estimated between 1 and 2 ft, 0.54 feet for o seems to&® be 
reasonable. So for the data used here, the measurements made 
by the Modified TDS Measurement System appear to be valid. 

3. Selecting Analyzer Settings 

Since the minimum time delay is 82 ms (Table 1), the 
maximum observation time top. which permits the HP3561A to 
resolve the direct and reflected signals is 80 ms. This 
corresponds to a 5-kKHz spans on the analyzer. As pointed 
out in Chapter II, choosing the maximum observation time 
gives the best frequency resolution on the analyzer. 

Figure 4.2 is the 3-D map output of Quick Analysis. 
The start frequency selected here is 500 Hz. Ideally, each 


5-kHz span in the sweep would be examined successively. 
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Figure 4.2: Quick Analysis 3-D Spectral Map 
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B. MEASUREMENTS 
The received Acoustic Level RL(f) at the face of the 


hydrophone is given by [Ref, l:p. 74]: 


RL(f) = ML(f)-CAL(f) -HS(f£)+IL(f)-47-AG 


where: 


ML(f) is measured level in (dB re 1V) 

CAL(f) is measured level of the calibration signal 
HS(f) is hydrophone sensitivity 

IL(f) is calibration insertion loss 


AG is the Auto Gain of the Amplifier. 


The levels ML(f) for these measurements are provided 
in outputs from the TDS _1 program in Appendix B. For 
selected frequencies used in this thesis, HS(f) and IL(f) 
are found in Table 2. These values are from the calibration 
records for the equipment provided by NUWES. The values for 
AG and CAL(f) used here are +40dGB and Od@B. Since KH sino is 
approximately unity at 1.2 kHz, some assumptions about 
Signal reflection at the surface can be made. For £f < gig 
kHz, the effect of specular reflection from the surface is 
determined by Equation (2-22). The upper and lower bounds 


of the acoustic pressure at the receiver are given by: 
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LTABEE 2 


HYDROPHONE SENSITIVITIES AND INSERTION LOSSES* 


Hydrophone 
Selected Upper Center Lower 
Frequency HS IL res eis HS is 
1 kHz alo G Oe O =—182.0 = 2 Oe Soc = 20m 
PeeKHZ eS dew Orne =eede > =2 0c elo. @ Oe 
Spek HZ aoe eee Ie Zak 2 Orme SO ie / One 
4 kHz -182.4 —20.0 =183 24 eee manG 2 2.3 =20.2 
S kHz eG. 3 = 2 Ole = Aceh ea 0, sue ee ene aa ze ONES 
*Data provided by NUWES 
XR 
Fmax = Pp(1 + xeR R.) (4-4) 
and 
XR 


Bom t ai |. 2ekhizZ othe wseffect of scattering and 
incoherent addition from the surface on the average acoustic 
intensity at the receiver is given by Equation (2-15) where: 


Me =e ee, leah (4-6) 


The ratio I./Iy is equal to the ratio of measured 


levels of both direct and reflected signal. Since the 
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Normal Analysis method uses data averaged over several 
Sweeps, these values Should pe Useq Ge ic a eameeie cet 
Multipath coefficients determined for the selected 


frequencies in Table 2 are displayed in Table 3. 


TABLE 3 


MULTIPATH COEPFICTEENTS AT SEEECTED FREGUE errs 


fl 2 type Coefficient 
1 Conenrehta( ley) DSS COs rae 
2 incoherent (My) ioe 
3 Incoherent we (May eee, 
4 incoherent (My) 1 eke, 
5 incoherent (M7) 1.18 


C. VARIABILITY 6: see seens 

Figure 4.3 showS a comparison of the results of 
Quick Analysis and Normal Analysis for the same conditions. 
Normal Analysis, as discussed earlier, makes an RMS average 
of 8 sweeps. In order to check the variability of results 
Quick Analysis was performed on six different sweeps. The 
energy averaged level of the measured band levels is given 


bys 


<X> = 10 10gi9 { }roLVL/10 yj ’ (4-7) 
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FIGURE 4.3: Quick Analysis and Normal Analysis Products 


(Output from TDS 1 Program) 
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where LVL is the measured level in a given frequency band 
and N is the number of measurements. Table 4 gives the mean 
level <X>, and "normalized" level of the standard deviation, 
10 log(S) - <X>, for several frequencies. 

The Rayleigh reflection parameter varies from about 
1 to 5.5 over the frequency span in the experiment. The 
energy averaged Quick Analysis measurements are close to the 
measurements obtained by the Normal Analysis method. This 
indicates that the Normal Analysis measurements provide a 
fairly good approximation to the average acoustic energy in 


each band of the measured signal. 


TABLE 4 


STATISTICAL COMPARISON OF SURFACE REFLECTED EE) Es 


(dB re 1V) 
Freq (HZ) Quick Analysis Normal Analysis 1O0log),4(2) > 22. 
<X> 
smerere) S26 = = Seo 
L500 = oto -14.6 =e 
2000 =o aie =k a 3 = Oecs 
206 =e en =2 0 el =r 
3000 = eee ele pets o 
4000 =2 470 ee - 4.0 
5000 -24.4 iD ee 
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V. DISCUSSION AND CONCLUSIONS 


A. BESecusSstOn 
1. Procedure 

Two factors which caused the greatest difficulty in 
analyzing data were the way in which the trigger was set and 
the choice of the number of sweeps for each data run. As 
discussed earlier, the Dynamic Signal Analyzer was triggered 
by signals recorded on the edge track of the data tape. The 
Honeywell 5600E magnetic tape recorder has a total of 14 
channels. This includes 13 data channels and one edge 
channel, typically used for voice comments. The edge 
channels can be recorded without disturbing the data 
recorded on the data channels. The edge tracks, located on 
the "edge" of the tape, sometimes suffer degradation caused 
by tape guides on the recorder. The trigger for the HP3314 
Arbitrary Function Generator was set to 4 seconds and 
recorded on the edge track. Several attempts were made to 
synchronize the trigger with the start of the LFM pulse. 
The data analyzed in this thesis required approximately 1.9 
second trigger delay. For the overlap processing used in 
the TDS program, only 80 percent of the LFM pulse was 
processed in Quick Analysis. Since the trigger was not 


synchronized to the start of the LFM pulse, 20 percent of 
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the signal between 8 kHz and 12 kHz could not be captured by 
the analyzer. 

The trigger level remained high for approximately 
1.7 seconds. This made Quick Analysis possible for parts of 
the sweep arriving during times of high trigger. 

Normal Analysis was not possible, however, because the 
analyzer continued to trigger and sample different portions 
of the LFM pulse. 

For further TDS measurements a sharp, short-duration 
trigger is necessary and should be recorded on tape with the 
data. The data Brekke used had a one second sweep period 
which was synchronous with the range time signals. A 1-Hz 
trigger was derived from the IRIG-B, time code generator 
translator. If a one second sweep period is maintained, the 
frequency span of the signal will have to be reduced to 
accommodate a slower sweep rate. If a slower sweep is used, 
such as the 4 second sweep used here, the advantage of a 
larger frequency span is gained. In this case, a trigger 
Signal synchronous with the start of the LFM pulse should be 
recouaed. 

The total data time for all three tests was 3.5 
minutes, with two minutes the maximum total time for slowest 
sweep rates. Normal Analysis requires 8 sweeps to make one 
measurement. For thirty sweeps, a maximum of 3 measurements 
can be made before the tape must be rewound, repositioned 


and started again. For the 5-kHz span used here 
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approximately 42 measurements are required and consequently 
the rewind process was repeated a minimum of 14 times. With 
five minutes of data recorded, the rewind process would be 
repeated only 5 times. Since the rewind process totals 
about 3 minutes the analysis time would be reduced 
Significantly. An additional benefit is that increased 
information available from the longer record could be used 
to determine the variability of the surface reflection. In 
future tests, the data recorded should be of a longer time 
at one sweep rate with 5 minutes of data accumulation as a 
minimum. 

2. Software 

The modifications made to the TDS 1 program have 
made it more flexible in making analysis of data. 
Eliminating a requirement for the plotter shown in Figure 
3.1 (b), would streamline the process considerably. 
Although HP Thinkjet graphic printouts are not as detailed 
or colorful as the plots of the HP7470A plotter, for 
perishable measurements the HP Thinkjet products are 
adequate. 

The Quick Analysis printouts given in Appendix B.1 
read in 4 ms increments; this corresponds to 16 ms of time 
in the Normal Analysis printouts, Appendix B.2. 

The program in current form can be used for routine 
measurements. Minor modifications need to be implemented 


for the changes just discussed, but personnel who would use 
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the products and conduct the measurements should determine 


the final form of the program. 


B;> -CONGLUS fGNS 
The Modified TDS Measurement System presented by Brekke 

(Ref. 1] provides useful information about an acoustic 
Signal propagating by direct and surface reflected paths. 
When the sea surface has some degree of roughness, 
measurements by this method permit the effect of the surface 
reflected signal and the direct signal to be considered 
separately. Further modifications to the measurement system 
are still required. In its current configuration, the 
system is cumbersome. The removal of the requirement for a 
plotter could streamline the analysis for routine tests. 

The HP3561A Dynamic Signal Analyzer can perform an 
energy average of a recorded bandwidth. This technique may 
be more accurate than the current method of recording only 
the peak value, particularly for Quick Analysis. In 
Normal Analysis the averaging over 8 sweeps, eliminates most 
of the variability so there should not be a significant 
difference between the average of peaks and the energy 
average. In Quick Analysis, however, the averaging over the 
bandwidth would improve the accuracy of measurements, 
particularly for the reflected signal at higher frequencies. 
The results obtained here are not conclusive. More data 
needs to be analyzed. The results reported here were made 


for a limited frequency range. Results for different sweep 
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rates and frequency spans should be investigated and 
compared. 

In any future tests, more extensive measurements and 
recording of environmental conditions need to be made. If 
reflected signals are distinguishable and if more than one 
is present, the operator may use TDS 1 to compare any of the 
reflected signals with the directly propagating signal. 
Future modification to the TDS program should allow for the 


analysis of more than one reflected signal. 


Bal 


APPENDIX A 


RANGE AND TEST DATA 


This list provides a guideline of information to be 


recorded at the range just prior to conducting tests using 


the Modified TDS Measurement System. 


lee 


ai 


4. 


ENVIROMENTAL 


(a) 
(b) 


(Cc) 


Wind-Speed and Direction 
Sound Velocity Profile 


Wave Heights 


CONFIGURATIONS CE TFs. 


(a) 
(b) 
(Cc) 
(d) 


EES 


Horizontal Ranges Between Source and Receivers 
Depths of Source and Receiver 
Range to Nearest Land 


Water Depth at Test Sight 


PARAMETERS 


(a) 
(b) 
(Cc) 


Calibration Signal 
Auto Gain Settings 
LFM Pulse 

- Frequency Span 

- Repetition Period 


- Number of Sweeps 


CALIBRATION RECORDS FOR RECEIVING EQUIPMENT 


(a) Hydrophone Sensitivities 


(b) Insertion Loss Data 
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APPENDIX B 
Tre ieee ASUREMENT RESULTS 


Meo: UPPER Rg DROPHONE 


ANALYIS METHOD: QUICK ANALYSIS 


DATE OP ANALYSIS: 12/1/86 


A SS SS hr a ce ec i ee ee ee ee 


RECORDED MAGNITUDE LEVELS (dB re 1V) 

DELAY(MS) DIRECT MAGN DIRECT PRO(HZ) REPL MAGN REPL PRO(HZ) NOISE MAGN 
1900 -9,33 500 19.75 oO -51.16 
1904 -9.21 563 19.75 oO -49.99 
1908 ~8.85 662 19.75 oO -49.69 
1912 -7.76 750 19.75 O -49.51 
1916 -9.28 875 19.75 oO -48.67 
1920 - 8.09 1013 19.75 oO -47.57 
1924 -~8.46 1125 19.75 oO -45.92 
1928 -6.39 1238 19.75 oO -43.42 
1932 -5.06 1300 -14.91 550 -42.00 
1936 -6.35 1363 -14.81 725 ~42.92 
1940 -9.85 1512 -13.85 787 -44.01 
1944 -11.53 1575 -14.13 eee -43.88 
1948 -9.96 1725 -12.80 1038 -42.29 
1952 -9.,93 - 1800 -10.42 1i25 -40.59 
1956 ~11.89 1913 -8.09 1263 -39.72 
1960 -12.27 2050 -7.14 1309 -40.14 
1964 -~11.24 2113 -9.77 1413 -41.08 
1968 -12.70 2263 -10.39 1438 -42.77 
1972 -11.76 2325 -17.93 1563 -43.72 
1976 -13.69 2438 -14.58 1788 -44.94 
1980 -15.24 2513 -11.66 1825 -44.79 
1984 -15.29 2675 -13.25 1862 -44.20 
1988 -13.84 2775 -15.29 2050 -43.45 
1992 -12.60 2875 -15.05 2150 -43.15 
1996 -11.92 2963 -13.69 2188 -43.06 
2000 -11.79 3075 -15.71 2325 -43.49 
2004 -11.76 3138 -21.23 2388 -43.11 
2008 -14.02 3250 -20.03 2525 -43.25 
2012 -14.14 3375 221,91 2662 -43.68 
2016 -13.78 3463 -22.64 2687 -43.72 
2020 -~ 14.54 3550 -19.82 2863 -44.14 
2024 —-16.51 3650 -20.56 2925 -44.03 
2028 ° - 16.61 3775 -20.28 3088 -44.07 
2032 - 16.57 3863 -19.28 3088 -43.96 
2036 -16.46 3950 -21.54 3175 -43.60 
2040 mi? , 92 4050 -22.95 3325 -43.76 
2044 -17.19 4188 -20.74 3513 -43.06 
2048 -16.93 , 4288 -19.11 3525 -42.96 
2052 17 07 4375 -25.74 3625 -42.80 
2056 -17.51 4475 ~18.46 3788 -43.04 
2060 - 17.06 4588 -19.91 3875 -43.03 
2064 - 17.76 4675 -21.03 3900 -43.03 
2068 -18.31 4763 -20.02 4113 -42.63 
2072 -18.96 4913 -~21.88 4150 -42.52 
2076 - 18.96 5012 -24.44 4250 -42.66 
2080 -18.49 5150 -25.28 4425 -42.89 
2084 -16.71 5262 -24.66 4463 -43.21 
2088 ~15.13 5313 -21.20 4625 -43.59 
2092 -~16.02 5388 -23.54 4662 -44.61 
2096 -16.85 5500 -25.07 4750 -44.73 
2100 -21. 43 5500 -27.18 4938 -44.39 
2104 19.75 0 -22.31 5038 -45.11 
2108 19.75 oO -21.15 5125 -45.06 
2999 14a =7¢6 a -20.63 5188 ~44.84 
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NPS: MIDDLE HYDROPHONE 


ANALYIS METHOD: QUICK ANALYSIS 


DATE OF ANALYSIS: 11/30/86 


EE ES Se ee hr 


RECORDED MAGNITUDE LEVELS (dB re 1V) 


DELAY(MS) OYRECT MAGN DIRECT PRO( HZ) REFL MAGN REFL FRO( HZ) NOISE MAGN 


——— — 


1930 -~6.60 563 19.75 O 47/81 
1934 7706 575 19.75 O -46.91 
1938 -4.90 750 19.75 O -46.20 
1942 -4,90 775 19.75 Oo -45.30 
1946 --6.5Z 925 19.75 O -43.57 
1950 -5.85 1075 19.75 O -41.11 
1954 = A662 1213 -10.46 537 =a 72 
1958 “23s 1275 =j1.74 613 -36.62 
1962 =e 1350 -12.51 775 BaG'.36 
1966 ~3.96 1400 -8.80 838 47-90 
1970 =7 469 1575 -10.68 912 -39,.35 
1974 maT S 1 i713 =o 04 1075 -—38.47 
1978 -7.56 1775 ST e22 1138 37.48 
1982 7.03 1250 Tes 1250 —37.80 
1986 ml.) 1950 ~8.40 1313 -39.74 
1990 - 10.60 2063 2S 1 1363 -41.76 
1994 - 12.05 2125 ~19, 36 1525 ~42.98 
1993 4.31 2225 -15.91 1713 -42.82 
2002 - 18.56 2375 214.30 1800 4252 
2006 ig 1G 2563 ea 77 1913 -41.81 
2010 mia.03 2625 -10.19 1975 -40.39 
2014 78 272 =A? . mA 2075 S40. 37 
2018 -9,89 2800 =e) 75 2162 ~41.09 
2022 9757 2888 = ae 2263 -42.49 
2026 Pet 22 2975 =26.56 2388 42 a7 
2030 62 3075 =o. 96 2513 4276 
2034 14.96 3250 -20.09 2600 240.19 
2038 =o 72 3325 26.20 2700 ~41.66 
2042 =e. 21 3438 73.33 2787 -40.95 
2046 2 787 3488 ~21.88 2838 -40.82 
2050 14.32 3625 -20.61 3075 -40.81 
2054 is. 76 3738 -16.74 3100 -~40.84 
2058 3,21 3825 -18.80 3238 -40.11 
2062 -, a4. 37 3938 =a 67 A559 39.70 
2066 - 13.54 4025 -16.05 3400 -39.76 
2070 12.61 4150 = 22.63 3475 -39.78 
2074 224 4713 Sirs 3.1. 3663 -39.95 
2078 - 11.85 4350 esas 3738 -38.92 
2082 -12.54 4400 ~15.81 3825 =e. 12 
2086 = 12 4 4550 =16.08 3938 age 79 
2090 e281 4625 gS 4088 -37.49 
2094 2780 4750 - 14.96 4188 37 41 
2098 e130 14 4850 Sa.) 4250 -~38.29 
2102 a7 76 4950 3215 4350 -38.96 
2106 -13.53 5038 =1'39 .67 4400 -~38.23 
2110 -12.81 5200 -15.56 4550 -38.56 
2114 2a 5225 =16.06 4588 -39.12 
2118 - 15.67 5288 -18.46 4763 -40.04 
2022 -15.49 5463 -19.72 4825 -40.48 
2126 =17.01 5500 - 22.65 5012 -41.58 
2130 19.75 0 20.71 5012 AQ a7 
2134 20161 563 22011 5200 -42.59 
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NPS: LOWER HYDROPHONE 


ANALYIS METHOO: QUICK ANALYSTS 


OATE OF ANALYSIS: 11/30/86 


i se ee 


RECORDEO MAGNITUDE LEVELS (dB re 1V) 


OELAY(MS) OIRECT MAGN OLIRECT PRO( HZ) REPL MAGN REFUL PRQO( HZ) NOISE MAGN 


a a rr RN CS ee 


1950 19.75 0 i9.75 0 oo co 
1954 he ae 725 19.75 0 =53.47 
1958 -11.40 800 19375 0 -52.60 
1962 -12.43 925 19.75 0 =o). 99 
1966 =12 62 1063 19.75 0 =30,, 91 
1970 =10.687 1150 Loe 72 0 -50.04 
1974 -8.40 1275 = 19), 97 537 -~48.96 
1978 aie 1325 =20.77 675 -48.13 
1982 — 9.92 1388 —19.33 800 -48.85 
1986 17 1500 Ole eed 850 ~49.74 
1990 -14.19 1688 72.6 1000 =49.13 
1994 =12.49 1738 -16.64 1075 -48.17 
1998 aaa 7 1813 -16.99 1150 -48.32 
2002 -15.20 1925 =19.58 i278 -49.24 
2006 eo +) 2050 -20.34 1363 ~49.33 
2010 -16.04 2138 -23.41 1475 -49.32 
2014 =17.56 2238 -21.40 1600 -49.20 
2018 ~18.68 2375 =19,.95 L725 -49.21 
2022 =19'. 70 2425 -16.20 eo -49.95 
2026 =o. 58 2600 7 47 1837 -49.89 
2030 =18.23 2700 =21,19 1950 -49.96 
2034 -16.74 2775 -26.30 2025 —@o.7¢ 
2038 Seo 7 2875 =2onee 2229 =49.12 
2042 15.73 2963 =27eol 2325 -49.24 
2046 =o, 32 3075 -30.01 2425 -49.31 
2050 = hei.o 7 3138 -26.34 2487 -48.66 
2054 =17 229 3300 -259.28 2588 ~48.14 
2058 16297 3400 -~24.01 2725 -47.54 
2062 =17.63 3463 =e ere 2850 =47.558 
2066 =197,,99 3600 -21.47 2950 -47.34 
2070 —19,01 3700 =l8251 3075 -46 .64 
2074 =19..a1 3800 -~16.43 3088 -46.81 
2078 1903 3900 =20,92 3188 -47 .17 
2082 -20.04 3975 -18.47 3325 -47.14 
2086 -20.06 4125 =21724 3362 —Ba 27 
2090 =19'.76 4225 =-26.53 3500 -47.24 
2094 =O ke 4325 —28 527 3663 -46.85 
2098 ~19.03 4438 -25.00 3763 -46 .63 
2102 =19701 ‘ 4525 =21.- 4) 3838 -46 .00 
2106 -18.64 4625 -20.66 3938 -45.49 
2110 =-19.25 4713 —23 eo 4088 -45.19 
2114 =. 9 14 1 4813 2 | 4150 35.55 
2118 - 19.48 4900 -22.86 4188 -46.10 
2122 “19,56 5025 =~29 527 4288 -46 .43 
2126 -18.28 5162 = 32.61 4500 -46.21 
2130 ho ae 5262 —20.05 4588 - 46.76 
2134 -17.76 5338 -28.89 4662 -47.82 
2138 -19.66 5400 ~29.68 4750 ~48.27 
2Vr42 =21.11 5500 —27 405 4913 -48.34 
2146 19275 Oo =22.19 5000 -48.96 
2150 19.75 0 —192e7 5075 -48.89 
2154 197 5 0 -19.65 5188 ~48.86 
2158 19.75 0 -20.04 5238 ~49.20 
2162 19.75 0 2a. §363 -49.70 
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ANALYTS METHOD: 


OATE OP ANALYSIS: 


DELAY( MS ) 


1900 
1916 
1932 
1948 
1964 
1980 
1996 
2012 
2028 
2044 
2060 
2076 
2092 
2108 
2124 
2140 
2156 
2172 
2188 
2204 
2220 
2236 
2252 
2268 
2284 
2300 
2316 
2332 
2348 
2364 
2380 
2396 
2412 
2428 
2444 
2460 
2476 
2492 
2508 
2524 
2540 
2556 
2572 
2588 
2604 
2620 
2636 
2652 
2668 
2684 
2700 
2716 
2732 


NPS: 


UPPER HYDROGPHONE 


NORMAL ANALYSIS 


RECORDED 


Loo 
19.75 
—9°. 14 
-10.98 
-13.94 
-11.66 
— soa 
19.75 
19.75 
19275 
= 776 1 
=10., 37 
19275 
-10.00 
~ 12,90 
193.75 
=12.03 
=13.90 
19279 
aia a 
-16.85 
-14.03 
15,76 
one. 
197,75 
19.75 
i375 
—15, 87 
oa el 
=o UG 
=<17 399 
19575 
=-16,95 
-18.65 
Loe75 
=. 
—lG@e27 
19.75 
—-17.26 
2h ane 
13275 
—17.02 
-18.99 
=17.62 
-18.86 
=—15 51 
-18.87 
=15 35 
=<17 757 
19.75 
-16.64 
19.75 
—-24.83 


12/1/86 





MAGNITUDE LEVELS 


DIRECT MAGN DOIRECT PROQ( HZ) 


775 
638 
787 
800 
963 


1388 
1738 


1738 
2050 


2088 
2350 


2325 
2613 
2438 
2787 
2975 


3275 
3500 
3362 
3650 


3688 
3950 


3975 
4263 


4288 
3925 


4600 
4888 
4688 
4963 
5262 
5012 
5288 
5500 


5475 


575 


56 


(dB re 1V) 


19.75 
19.75 
19.75 
19.75 
19-72 
19:75 
i. 75 
19.75 
19275 
19.75 
~laoo4 
=-15.72 
19.75 
=—15°.14 
=12,39 
19.75 
—9.50 
71367 
19.75 
=14;15 
-13.45 
-14.45 
-11.38 
=1957 9 
19.75 
19.75 
59575 
—22.67 
-20,;52 
-20.24 
=17 -O2 
19.75 
-~20.10 
-18.58 
19.75 
-18.886 
-29.36 
19.75 
-23.68 
=25,22 
19.75 
=2 17.90 
~19.34 
=—21,09 
-20.25 
~—eeno 
—-24.54 
-22.95 
=22.17 
19.75 
-24.00 
-22.76 
=2 Si) 2 


REPL MAGN REPL PRQ( HZ) 


ooo oo 00 000 


725 
1025 


1063 
i325 


i375 
i7i3 


1688 
1913 
1775 
975 
2225 


2588 
2825 
2662 
2988 


3000 
3225 


3250 
3513 


3650 
3925 


3900 
4188 
4012 
4238 
4538 
4288 
4563 
4850 


4725 
5062 
5225 


NOISE MAGN 


-47.29 
~49.07 
-46 .88 
-47.79 
=A aaa, 
-46.11 
- 46.98 
—-46.98 
~46.98 
-50.59 
=o oe 
-40.88 
-40.88 
-40.72 
-~38.68 
-38.68 
=—2o oS 
=41.26 
-41.26 
-41.00 
=39. 61 
-41.69 
-39.08 
Ag. 307 
ig PARE 
—42.37 
Sg PSOE), 
~40.30 
—=42.85 
tm ea 
=—39.05 
=—29 205 
=41.72 
-38.69 


-41 ° 


| 
SReTe 


-40.76 
-40.76 
Io. oT 
=38).56 
=39 92 
=39:.87 
~40.13 
-40.12 
-40.37 
—-41.76 
-41.76 
-41.09 
-40.32 
-~43.41 


ANALYIS 


DELAY( MS ) 


1930 
1946 
1962 
1978 
1994 
2010 
2026 
7042 
2058 
2074 
2090 
2106 
2122 
2138 
2154 
2170 
2186 
2202 
2218 
2234 
2250 
2266 
2282 
2298 
2314 


ZIAWM 


tw 


2346 

2362 

2378 

2394 
2410 

2426 

2-42 

2458 

2474 
2490 

2506 
2522 
2538 
255% 
2570 
2586 
2602 
2618 
2634 
2650 
2666 
2682 
2698 
2714 
2730 
2746 
2762 


METHOD: 


DATE Of ANALYSTS: 


NE So 


NORMAL ANALYSIS 


RECORDED MAGNITUDE LEVELS 


DIRECT MAGN 


7 oe 

-5.68 

-7.43 

-6.07 

-6.92 

- 2.66 

SST 

-4.24 
—11 41 

-8.44 

ats od 1 

i. 75 
Bele 2). 6 1 
~15.45 
=12.28 
- 16.65 
meee S 
-~20.14 
-14.61 
- 16.73 
aes Jo 
gaa. 70 
a1. 45 
=14:29 
= 14.28 
—13',08 
ad 65 
= 157, 95 
12.60 
seo. 0 2 
~18.46 
3°, 390 
> LS eto 
14, Bl 
—lf.e2 
— 13.439 
=o. Lo 
-13.6)1 
- 14.54 
712.97 
-14.72 
- 15.36 
- 13.83 
ot. 22 
- 14.23 
—13.65 
- 17.83 
= 15,08 
Ge 7k 
- 47.10 
- 46.87 
=45,35 
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MIDDoUE HYDRPOPHONE 


12/1/86 


DIRECT PROQ( HZ) 


500 
750 
963 
787 
1100 
1313 


1388 
1713 
1400 
1775 


1862 
2150 
1988 
2238 
2138 
2450 
2737 
2550 
2800 
2963 
2787 
3050 
3325 
3050 
was 7 
3563 
3387 
3738 
3925 
3775 
4037 
3850 
4150 
4425 
4238 
4475 
A725 
4538 
4813 
5062 
4850 
5188 
5012 
5238 
5500 
5262 
5500 
3300 
4163 
3288 


Bi 


(dB re 1V) 


REPL MAGN REPL PROQ( HZ) 


Lo. 75 

Lo 

13 79 

3.75 

29.75 
-13.44 

T3279 
=1).57 

ar at 
-12.00 
~10.34 

3.75 
-10.84 
-14.66 
-7.44 
=15. 10 
aes 93 
een ok 
23. 90 
-~14.09 
ooo. 
-30.08 
=—18). 14 
—30.94 
-24.50 
-26.08 
eine 
= 217.00 
-22.94 
eis Ld 
-18.28 
- 21.60 
-~17.45 
-16.75 
-~20.62 
—17 95 
~16.3% 
-21.49 
-19.42 
leo 
- 19.53 
- 20.28 
~17.80 
—lveee 
a sO 
~18.16 
-22.87 
=18.82 
-22.69 
ee 7S 
~23.76 
-24.53 
SO.77 


OoO0o00 0 


750 


825 
1275 
S12 
1213 


1313 
1575 
1375 
1688 
1475 
1775 
2150 
1925 
2113 
2400 
2088 
2463 
2725 
2425 
2712 
3013 
2800 
3125 
3463 
gi13 
3387 
e312 
3538 
3825 
3650 
3950 
4150 
3925 
4250 
4413 
4300 
4S50 
4388 
4650 
4975 
4700 
5012 
§212 
5138 
5350 
5188 


-43.05 
-40.55 
—37, 96 
-42.38 
-37.34 
-34.% 
~34.% 
— 39.07 
—39.05 
-35.40 
-33.98 
=—33.98 
—-36.46 
=—I9 702 
—35.. 96 
=4aG. 95 
=I co 
=—36.94 
-38.87 
=35 19 
=a» 70 
-40.76 
—39.795 
—-34.28 
-39.68 
-40.02 
~37.96 
=37 762 
—39.30 
-36.58 
-~36.68 
~36.10 
-36.16 
—38.,12 
-~36.54 
- 36.87 
7 Od 
=3505 
-37.45 
36 ..7 1 
-34.68 
-38.15 
-36.48 
~34.21 
-36. 44 
-35.15 
~39.41 
-37.40 
-36.37 
-40.03 
-40.44 
-37.55 
-39.64 


NOISE MAGN 


NPST: LOWER HYDROPHONE 


ANALYIS METHOD: NORMAL ANALYSIS 


DATE OP ANALYSIS: 11/30/86 


RECORDED MAGNITUDE LEVELS (dB re 1V) 


DELAY(MS) OILRECT MAGN OIRECT FPRQ( HZ) REPL MAGN REPL PRO( HZ) NOISE MAGN 


1950 -14.98 662 19.75 oO -51.50 
1966 ~14.68 912 19.75 oO -48.51 
1982 -11.24 1188 19.75 oO -48.03 
1998 -13.12 975 19.75 oO -50.42 
2014 ~9.93 1263 -21.03 575 -47.79 
2030 ~14.56 1388 -18.47 813 -47.43 
2046 -9.72 1825 -21.52 675 -47.14 
2062 -14.58 1563 -16.53 875 -46.68 
2078 -18.08 1775 -19.52 1250 ~44.77 
2094 ~14.41 1700 -—20.90 1050 -46.75 
2110 -15.23 1900 -15.40 1313 -45.23 
2126 19.75 oO 19.75 0 -45.23 
2142 -16.86 2075 -13.21 1388 -44.71 
2158 19.75 oO 19.75 0 -44.71 
2174 -17.68 2213 -17.51 1475 -46.16 
2190 -21.85 2525 -15.29 1813 -46..43 
2206 -18.56 2288 -22.28 1700 -44.55 
2222 ~19.50 2625 -20.51 1913 -46.75 
2238 19.75 oO 19.75 oO -44.94 
2254 -17.98 2725 -20.10 1975 -46.54 
2270 -16.99 2963 -24.49 2263 -47.05 
2286 ~16.51 2825 -23.31 2075 -46.45 
2302 ~16.05 3125 -~25.68 2475 -46.14 
2318 --16.32 3050 -27.78 2400 -46.23 
2334 -17.49 3350 -26.05 2712 -45.83 
2350 -18.22 3250 -29.23 2613 ~45.74 
2366 -23.30 3588 -23.90 3025 -46.55 
2382 -18.01 3450 -26.33 2813 -45.79 
2398 -20.50 3763 -18.76 3075 -45.33 
2414 -20.27 3650 -22.28 3013 -45.67 
2430 -21.14 3963 -19.63 3225 -44.09 
2446 -20.07 3838 -21.37 3113 -45.34 
2462 -20. 42 4188 -27.98 3438 ~44.97 
2478 —21.12 4063 -23.38 3350 -45.16 
2494 ~21.38 4388 -28.26 3788 -44.19 
2510 19.75 oO 19.75 oO -44.19 
2526 -~20.11 4488 -25.96 3825 -44.24 
2542 ie. 4775 -25.18 4050 -42.58 
2558 - 19.76 4575 -25.24 3888 -43.69 
2574 - 19.69 4900 -25.83 4200 -43.90 
2590 - 20.20 4750 -26.31 4125 —-4%.22 
2606 ~18.95 5100 -27.98 4388 ooo 
2622 -20.06 5000 -24.44 4338 -42.88 
2638 -16.85 5288 -24.01 4650 - 44.05 
2654 -17.63 5225 -25.84 ae hee: 
2670 -19.39 5438 -24.14 4775 -44.95 
2686 -20.82 5413 -26.16 4775 -46.11 
2702 19.75 0 =27 42 5025 -47.59 
2718 19.75 Oo -29.21 5188 -47.33 
2734 19.75 Oo 19.75 oO -47.33 
2750 19.75 oO —25.72 5188 -46.33 
2766 19.75 oO =27.77 5500 -47.66 
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